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1 Introduction 



Consider the one dimensional liquid-gas two-phase model with viscosity which can be 
written in Eulerian coordinates as (cf. OH]): 



dt[a g pg] + d x [a g p g Ug} = 0, 
d t [atipi] + d x [aipiui] = 0, 
k dt[a g p g u g + aipiui] + d x [a g p g u 2 g + aipiuf + P] = -q + d x [ed x u r , 



(1.1) 



OigUg + a\Ui and the unknown variables ce g ,cti £ [0, 1] denote volume fractions 



where u m i x — ^g^g 
satisfying the fundamental relation: 



«„ + «! = 1. 



Cl-2) 



Furthermore, the other unknown variables p g , pi, u g , ui denote gas density, liquid density, 
velocities of gas and liquid respectively, whereas P is the common pressure for both phases, 
q presents external forces, like gravity and friction, and e > denotes viscosity. 

We focus on a simplified model as in O H] obtained by assuming that fluid velocities 
are equal, i.e., u g = ui = u and neglecting the external forces, i.e., q = 0. In addition, we 
neglect the gas phase effects in the mixture momentum conservation equation (1.1)3. This 
motivation is from the fact that liquid density is much higher than the gas density, generality 
speaking, pi/p g = O(10 3 ). Thus, we can obtain the following simplified model: 

d t [a g pg] + d x [a g pgu] = 0, 

dt[aipi] + d x [aipiu] = 0, (1.3) 

k d t [aipiu] + d x [aipiu 2 + P] = d x [ed x u]. 

As in [3], we assume the liquid is incompressible, i.e., pi = constant and the gas is poly- 
tropic, i.e., 

P = Cp], 7 >1, 0. (1.4) 

Let 

n = a g p g , m = aipi. (1.5) 
By (1.2), (1.3), (1.4) and (1.5), we have 

d t n + d x (nu) = 0, 

dtm + d x {mu) = 0, (1.6) 
k d t {mu) + d x (mu 2 + P) = d x (sd x u), 



and 



P(n, m) 



CP, 



n 



pi-m 



A 



n 



pi-m 



where A = Cpl. Moreover the viscosity coefficient is taken as the form (cf. [3 S]): 



e = e(n, m) = B 



rv 



(pi - m)? 



B > 0, p > 0, 



(1.7) 



2 



or 



e = e(m)=B- —-, B > 0, /3 > 0. (1.9) 

(pi - m)P+ 1 

In the following, without loss of generality, we consider only the form (1.8). 

In this paper, we will consider the following two classes of the free boundary value prob- 
lems of the system (1.6): 

(1) The initial masses connect to vacuum discontinuously: 

The boundary conditions are given as 



"" U! ^u(a(t),t), t > 0, 




and 

db(t) 



u(b(t),t), t>0, 



(1.10) 



(— P(m, n) + e(n,m)d x u)(a(t) + ,t) = 0, 
(-P(m, n) + e(n, m)d x u)(b{t)~ ,t) = 0, t > 0, 
and the initial data are given as 

n(x, 0) = rio(x) > 0, m(x, 0) = mo(a?) > 0, u(x,0) = no(x), x G [a, b]. (1-H) 

(2) The initial masses connect to vacuum continuously: 
The boundary conditions are given as 

n(a(t),t) = n(b{t),t) = 0, m(a(t),t) = m(b(t),t) = 0, t > 0, (1.12) 

and the initial data are given as 

n(x, 0) = no(x) > 0, m(x,0) = mo(i) > 0, u(x, 0) = uo(x), x G (a, 6), (1-13) 

and no(a) = rio(6) = mo(a) = mo(6) = 0. 

Here — oo < a < b < oo, a(i) and 6(t) are the free boundaries defined by 



;i.i4) 



dt v w ' y ' (1.15) 
6(0) = 6, 

which are the interfaces separating the gas from the vacuum. 

Let's first review some of the previous works in this direction. When the viscosity coef- 
ficient e was taken as the form (1.9) and the initial masses connected to vacuum discontinu- 
ously, Evje and Karlsen in [3] got the global existence and uniqueness of weak solutions when 
ft G (0, |) by energy method. This result was later generalized to the case when (5 G (0, 1] by 
Yao and Zhu in [17j . When the viscosity coefficient e was taken as the form (1.8) and the ini- 
tial masses connected to vacuum continuously, Evje, Flatten and Friis in [1] proved the global 
existence of the weak solutions under some approximate assumptions on no(x), mo(x), ^0(2;), 

and = co(x) when f3 G (0, A). This result was later generalized to the case /3 G (0, 1) 

too(x) 

by Yao and Zhu in [16] when the initial masses connected to vacuum continuously and the 
viscosity coefficient was a constant. 
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It is necessary for us to illustrate that the main methods used to obtain our results are 
similar to those in [23\ [7J, because we have used the variable transformations as in [31 E], 
by which we can rewrite our problem into (2.7)-(2.12) similar to the model in single-phase 
Navier-Stokes equations. In view of this, let's review some of the relevant works about 
single-phase Navier-Stokes equations with density-dependent viscosity and vacuum. When 
the initial density connected to vacuum discontinuously, the global existence of weak solutions 
for isentropic flow was obtained by Okada, Matusii-Necasova, Makino in [15] for p,(p) = p e , 
< 9 < 3, by Yang, Yao and Zhu in |19j for < 9 < \ an d by Jiang, Xin and Zhang 
in [H] for < 9 < 1. Qin, Yao and Zhao in fTS] extended the results in [TS1 H2J [H] to the 
case < 9 < 1. Recently, Zhu in [23J investigated the asymptotic behavior and decay rate 
estimates about the density function p(x, t) by overcoming some new difficulties which came 
from the appearance of the boundary layers. When the initial density connected to vacuum 
continuously, the local existence of weak solutions was obtained in |20j by Yang and Zhao. 
The global existence of weak solutions was given in [21] by Yang and Zhu for < 6 < |, and 
later was improved in [22] for < 9 < ^ by Vong, Yang and Zhu, in [5] for < 9 < | and 
in [6] for < 9 < 1 by Fang and Zhang. Guo and Zhu in [7J firstly studied the asymptotic 
behavior and the decay rate of the density function p(x, t) with respect to the time t for any 
9 > based on the following new mathematical entropy inequality, which was obtained first 
by Kanel in [9] for one-dimensional case and Bresch, Desjardins, Lin and Mellet, Vasseur for 
multi-dimensional case, cf. [T| |2"1 IT2]: 



where C is a uniform constant independent of t. 

In this paper, we show that the masses n and m tend to zero as time goes to infinity. 
Moreover, we can obtain a stabilization rate estimates of the mass functions for any /3 > 
as t — > 00. 

The rest of this paper is organized as follows. In Section 2, we reformulate the two free 
boundary value problems (1.6), (1.10), (1.11) and (1.6), (1-12), (1.13) into the two fixed 
boundary value problems by introducing the Lagrangian coordinates and using the variable 
transformations. Then we state the main theorems of this paper. In Section 3, we derive some 
crucial uniform estimates for studying the asymptotic behavior and the decay rate estimates 
about the mass functions. In Section 4, the decay rate estimates on the mass functions will 
be given by introducing a new function w(x, t) in |13] by Nagasawa. 

2 Reformulation of the problems and the main results 

To solve the two free boundary problems above, it is convenient to convert the free 
boundaries to the fixed boundaries by using Lagrangian coordinates. To do this, let 




} 



dx 




T = t. 
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Then the free boundaries x = a(t) and x = b(t) become £ = and £ = f^jty m(y, t)dy = 
J m o{y)dy by the conservation of mass, where mo(y)dy is the total liquid mass initially. 
We normalize mo(y)dy to 1. 

Hence in the Lagrangian coordinates, the two free boundary problems (1.6), (1.10), (1.11) 
and (1.6), (1.12), (1.13) become 

n T + (nm)u£ = 0, 

m T + m 2 U£ = 0, (2.1) 
u T + (P(n,m))^ = (e(n,m)mu^)^, 

with the boundary conditions (corresponding to the initial masses connect to vacuum discon- 
tinuously) 

P(n,m) = E(n,m)v,(, at £ = 0,1, r > 0, (2.2) 

and the initial data 

n(£,0)=n o (£)>0, m(£, 0) = m (£) > 0, u(£, 0) = u (£), £€[0,1], (2.3) 

or with the boundary conditions (corresponding to the initial masses connect to vacuum 
continuously) 

n(0, r) = n(l, r) = 0, m(0, r) = m(l, r) = 0, r > 0, (2.4) 
and the initial data 

n(£,0) = n o (£)>0, m(£, 0) = m (£) > 0, u(£, 0) = u (£), £€(0,1), (2.5) 

and no(0) = rao(l) = ^o(0) = w-o(l) = 0. 
Here 

, e(n,m)m = - -5— r = EYn, m), /3 > 0. (2.6) 

pi-mj {pi-m)P +l 

Here we have assumed ^4 = B = 1 in (2.6) for simplicity. 

In the following, we replace the coordinates (£, r) by (x,t). Introduce the variables (cf. 

c=— , Q(m) = = — ^>0. 

m pi — m 1 — ai 

Form the first two equations of (2.1), we get 

nt n mnu x nm 2 

c t = ~m t = 1 ttU x = 0, 

m m z m m z 

and 

/ m \ ( 1 m \ 

Qi m )t = ( - — — ) = ( z — ~ + tz — ZZ^2 m * 



pi-mj t \pi-m (pi-mY 



Pi Pim 2 

m t = -7 To u i 



(pi - m) 2 (pi - m) 2 

-piQ(m) 2 u x . 
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Then we can rewrite the initial boundary problems (2.1), (2.2), (2.3) and (2.1), (2.4), (2.5) 
into the following forms: 

( d t c = 0, 

d t Q(m) + Pl Q{m) 2 d x u = 0, (2.7) 

k d t u + d x (P(c,m)) = d x (E(c,m)d x u), 

with the boundary conditions (corresponding to the initial masses connect to vacuum discon- 
tinuously) 

P(c,m) = E(c,m)u x , at x = 0, 1, t > 0, (2.8) 

and the initial data 

c(x,0) = cq(x) > 0, Q(m)(x,0) = Q(m ){x) > 0, u(x , 0) = u {x) , s€[0,l], (2.9) 

or with the boundary conditions (corresponding to the initial masses connect to vacuum 
continuously) 

c(0, t) = c(l, t) = 0, Q(m)(0, i) = Q(m)(l, t) = 0, t > 0, (2.10) 
and the initial data 

c(x, 0) = c (x) > 0, Q(m)(x, 0) = Q(m )(x) > 0, u(x, 0) = u (x), x G (0, 1), (2.11) 

and c (0) = c (l) = Q(mo)(0) = Q(m )(l) = 0. 
Here 

P(c,m) = f 71 V = c 7 Q(m) 7 , E(c,m) = me(n,m) = c^Q(mf +1 , (3 > 0. (2.12) 
VP/ ~ m / 

Throughout this paper, our assumptions on the initial data and /3, 7 are stated as follows: 

(Ai) inf no(x) > 0, sup no(x) < 00, inf mo(x) > 0, sup mo(x) < pi; 
xe[o,i\ xe[o,i] xe[o,i] ie[o,i] 

(Ai)' There are positive constants K\, K2, K$ and such that K\4>(x)~ < mo(x) < 
K 2 (j)(x)^ < pi, K 3 <p(x) a < n (x) < K 4 <f>(x) a , where <p(x) = x(l - x), < a < 1. In 
particular, this implies that there exist positive constants Ci, C2, shuch that C\(f){x)^ < 
= S$ < G.*(*)t; 

(A2) no(x) € L 2ra ([0, 1]) for any given positive integer n satisfying n > 2 

(A 3 ) ((coQ(m ))^GL 2 ([0,l]); 

(A4) /3>0, 7 >l + /3. 

Now we give the following definition of weak solution: 

Definition 2.1 (Weak solution). We call (n(x,t),m(x,t),u(x,t)) a global weak solution 
to the initial boundary value problems (2.1), (2.2), (2.3) or (2.1), (2.4) (2.5), if the following 
estimates hold for any t > 0, 

n, m, u G L°°([0, 1] x [0, +00)) H C 1 ([0, +00); iJ 1 ([0, 1])), 
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E(n,m)u x £ L°°([0,1] x [0, +00)) n ([0, +00); L 2 ([0, 1])), 
< n(x,t) < pi sup Co, 

X6[0,l] 

and 

< m(x, t) < pi. 
Furthermore, the following equations hold: 

rit + mnu x = 0, mt + m 2 u x = 0, 

(n, m)(x, 0) = (no(x), mo (a;)), for a.e. x G [0, 1] and any t > 0, 

and 

/ / (utpt + (P(n,m) — E(n,m)u x )ip x )dxdt + / uo(x)9?(x, 0)dx = 0, 
jo Jo Jo 

for any test functions 93 G C :x:> (O) with Q = {(x, i) : < x < 1, t > 0}. 



In what follows, we always use C (and C n ) to denote a generic positive constant depending 
only on the initial data (and the given positive integer n), but independent of t. 

We now state the main theorems in this paper as follows: 

Theorem 2.2 (The asymptotic behavior of the mass functions). Under the assump- 
tions (Ai), (A2), (A3) and (A4) (corresponding to the boundary conditions (2.8)) or (Ai)' , 
(A2), (A3) and (A4) (corresponding to the boundary conditions (2.10)), let (n(x,t), m(x,t), 
u(x,t)) be a global weak solution to the initial boundary value problem (2.1), (2.2), (2.3) 
or (2.1), (2.4) (2.5). Then we have the following asymptotic behavior of the mass functions 
n(x,t), m(x,t) 

lim sup n(x,t) = 0, (2.13) 
*^°°xe[o,i] 

lim sup m(x,t) = 0. (2-14) 

*^°°x6[0,l] 

Furthermore, we can get the decay rate estimates of the mass functions n(x,t), m(x,t) 
as follows: 

Theorem 2.3 (The decay rate of the mass functions). Under the assumptions of 
Theorem 2.2, let (n(x,t), m(x,t), u(x,t)) be a global weak solution to the initial boundary 
value problem (2.1), (2.2), (2.3) or (2.1), (2.4) (2.5). Then the following decay rate estimates 
on the mass functions n(x,t), m(x,t) hold: 

(i) Under the boundary conditions (2.8), if < f3 < 1 or (3 > 1, ^ > 2, then 

n(x,t), m(x,t) <C(l + t)~^+^, (2.15) 



for any x G [0, 1] 



If P = 1 or > 1, 3=i < 2, we have 



n(x,t), m(x,t) <C(l+t) ^w(ln(l + t))^W, (2.16) 
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for any x G [0, 1]. 

(ii) Under the boundary conditions (2.10), if < j3 < 1 or /3 > 1, > 2, then 

n(x,t), m(x,t) < C(l + t)~y-°+w , (2.17) 

for any x € [0, 1]. 

\ 1 

7-/3 



If = 1 or p > 1, ^ < 2, we have 



n(x,t), m(x,t) < C(l + t) n-i+*P(]n(l + t))n-^P, (2.18) 
for any x 6 [0, 1], where # is defined by Lemma 4.1. 

3 A priori estimates and the asymptotic behavior of the mass 
functions 

In this section, we will give some useful uniform a priori estimates of the solutions with 
respect to the time t. Then we study the asymptotic behavior of the mass functions n(x, t) 
and m(x, t) by using these uniform a priori estimates. 

3.1 Uniform a priori estimates 
Lemma 3.1 (Some identities). 

c(x,t) = c (x), (3.1) 
d f x d f 1 

di Jo u ( y,t ) dy = ~dt J u ( y ^ dy ' ^ 3 - 2 ) 
(cPQ(mf +1 u x )(x,t) = (ciQ(my)(x,t) + f u t (y,t)dy = (JQ(m,y)(x,t) - C u t (y,t)dy, 

JO Jx 

(3-3) 

-L( c PQ(mf)(x,t)+ [ c<Q(my(x,s)ds = -?-4Q(mo)P - f f u t (y,s)dyds. (3.4) 
PPi Jo PPi Jo Jo 

Proof. These identities can be obtained directly from (2.7). 

Lemma 3.2 (Basic energy estimate). Under the conditions in Theorem 2.2, the following 
energy estimate holds: 

u 2 + 



2 Pib 

1/1 „7 



—Q(my l {x,t)]dx+f [ c fj Q(mf +1 u 2 x dxds 
1) / Jo Jo 



-I (^ + M^T)^-')' fa 

< C. (3.5) 

Proof. Multiplying the second and the third equations of (2.7) by c 1 Q{m) 1 ~ 2 and u, and 
integrating the resulting equations with respect to x over [0,1], we get 

Sr(T + M^T) <3( " !)7 i' il ' + <, ' P(c - m>) 

1 f 1 

= (E(c,m)u x u) — / E(c,m)u 2 dx. (3-6) 
Jo 
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Using the boundary conditions (2.8) or (2.10) we get 

j t £ (y + ^^y QM 7 " 1 ) ^ + cPQ(mf +1 u 2 x dx = 0. (3.7) 



Then integrating it with respect to t over [0, t], we get (3.5). 
The proof of Lemma 3.2 is completed. 

Lemma 3.3 (The Uniform upper bound for the cQ(m)). Under the conditions of 
Theorem 2.2, we have for any xe [0,1], t > 0, 

< (cQ(m))(x,t) < C. (3.8) 

Proof. By (3.4), we have 

-L(cPQ(mf)(x,t)+ f ciQ(my(x,s)ds = J-^Q(m f - f \u(y,t) - u (y))dy 
PPi Jo PPi Jo 



<C+ [ \u(y,t)\dy+ [ \u {y)\dy 
Jo Jo 

< c + u2 dy) 2 + (J q 4dy 



< c. 

Here we have used Lemma 3.2 and the assumptions (^4i), (Ai)' , {A2). 
The proof of Lemma 3.3 is completed. 

Corollary 3.4. For any x <G [0, 1] and t > 0, 

/ c'Q{my(x,s)ds < C. (3.9) 
Jo 

Lemma 3.5. Under the boundary conditions (2.8), we have for any t > 0, 

Q(m)(d,t) = Q(m )(d)(- i~"g~~~~ : „ -) \ d = Q, 1. (3.10) 

Proof. By (3.2), we have 

d f 1 

jtj u ^ d y = °- 



By taking x = 1 or x = in (3.4), we have 

ft 

Pp, 

Since 7 > j3, the integral equation (3.11) yields 



1 -4(d)Q(mf(d,t)+ [\m)Q{my(d,s)ds = -l-4{d)Q(mo) l3 (d), d = 0, 1. (3.11) 
W Jo PPz 



-cg^QH^^^Hfd,^ + c2(<l)QK((l,() = 0, (3.12) 
Pi 

which implies (3.10) by solving the ordinary differential equation (3.12). 



The proof of Lemma 3.5 is completed. 

Corollary 3.6. There exist positive constants C\ and C2 such that for any t > 0, 

Ci(l + t)"^ < Q(m)(d,t) < C 2 (l+t)~^. 

Lemma 3.7. For any positive integer n in (A2), we have for any t > 

f u 2n dx + n{2n-l) [ [ cf ] Q{mf +l u 2n - 2 u 2 x dxds < C n , (3.13) 
Jo Jo Jo 

where C n is a positive constant depending on n, but independent of t. 

Proof. Multiplying the third equation of (2.7) by 2nu 2n_1 and integrating the resulting 
equation with respect to x over [0, 1], we have 

j t ^ u 2n dx + 2n{u 2n - l P{c,m))^ - 2n(2ra - 1) ^ c 1 Qim) 1 u 2n ~ 2 u x dx 

1 /•! 

= 2n(n 2n - 1 J E(c,m)n a; ) - 2n(2n - 1) / c^m^+V™" 2 ^^. (3.14) 

io 

Using the boundary conditions (2.8) or (2.10), we have 

^ C u 2n dx + 2n(2n - 1) C ^Q{mf +1 u 2n - 2 u 2 x dx = 2n(2n - 1) C c~ ( Q(mf< u 2n ~ 2 u x dx. 
dt Jo Jo Jo 

(3.15) 

Integrating (3.15) with respect to t over [0,t], we get 

/ u 2 "(ix + 2n(2n - 1) / / c^Q(m)^ + V n -\ 2 dxds 
jo Jo Jo 

= [ u 2 n dx + 2n(2n-l) [ [ c- 'Q(m)" 'u 2n ~' l u x dxds . (3.16) 
Jo Jo Jo 

Applying Cauchy-Schwarz inequality to the last term in (3.16) yields 

/ u 2n dx + n(2n - 1) / / ^Q{mf +1 u 2n - 2 u 2 x dxds 
J \ J \ J \ (3.17) 

< ( ul n dx + n{2n-l) I I c 2 ^^ Qimf^^u^dxds. 
Jo Jo Jo 

Now we estimate the last term on the right-hand side in (3.17) as follows: 



n(2n-l) [ [ c 2 ^Q(m) 2 ^- x u 2n - 2 dxds 
Jo Jo 

= n(2n-l) f f cl + ^-PQ(m)% + ^-P- 1 c^Q(m) I ^r'- / u 2n - 2 dxds 
Jo Jo 

<{2n-l) [ [ c 1+n ^- p) Q{my +n ^- p - l) dxds + {n-l){2n-l) [ [ c 7 Q(m) ~<u 2n dxds 
Jo Jo Jo Jo 

= (2n - 1) / / (cQ(m)) n ^-P-Vc n (cQ(m)y'dxds + (n - l)(2n - 1) / / c / Q(m^u 2n dxds 
Jo Jo Jo Jo 

<C I max(cQ(m)) 7 ds + C [ max(cQ(m)) 7 / u 2n dxds 
Jo [o>i] Jo [0-1] Jo 

<C + C f max(cQ(m)) 7 [ u 2n dxds. 
Jo I - 1 ] Jo 



(3.18) 
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Here we have used the Young inequality ab < ^ + w h ere | + \ = 1> Pi Q. > 1> a ; ^ > 0, 
the assumption (A4), (3.8) and (3.9). 

Substituting (3.18) into (3.17), we have 

u 2n dx + n(2n - 1) f C ] Q{mf +1 u 2n ~ 2 u 2 x dxds 
Jo Jo Jo 

< C + C I max(cQ(m)) 7 / u 2n dxds. (3.19) 
Jo Ml 7o 

By (3.19), (3.9) and GronwalPs inequality, we have 

/ u 2n dx < C n , (3.20) 
Jo 

where C n is a positive constant depending on n, but independent of t. 

(3.19) and (3.20) show that (3.13) holds and this completes the proof of Lemma 3.7. 

Lemma 3.8. We have the following uniform estimate on the derivative of the function 

cQ(m), 

f ({cQ{m)f) 2 dx+ [ [ ((cQ(m))^) 2 dsdx <C. (3.21) 
^ ' x Jo Jo ^ ' x 



Proof. From (2.7), we have 

{{cQ{m)f) xt = {P{cQ{m)f-\cQ[m)) t ) x 
= (jScPQimf^Qim)^ 
= -ppi{c^Q{mf +1 u x ) x 

= -l3 Pl {u t + P{c,m) x ). (3.22) 

Multiplying (3.22) by {{cQ{m)Y) x , integrating the resulting equation over [0, 1] x [0, t] and 
integrating by parts, we get 

- y ((cQMy 3 )^^- 1 ((cogN)^^-^, y ^(cQH^dx 

+Ppi J u ((c Q(m ))^ ^dx + f3pi J J u{{cQ{m)fy^dsdx 

(3.23) 
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Substituting (3.22) into (3.23), we have 



+P Pi / «o ( (coQimo)) 13 ) dx-(/3pi) 2 / / uu t dsdx 
Jo v /;c jo io 

-(/3 Pi ) 2 / / u{{cQ{m)V) x dsdx 
Jo Jo 

= \J ({c Q Q{m )Yyjx-i3p l J u[{cQ{m)f)jx 
+P Pl f n ((c Q(m )f) Jx f u 2 dx + ^ 



\2 



-(/3^) 2 / {(cQ(m)y(l,s)u(l,s)-(cQ(m)y(0,s)u(0,s)}ds 
Jo 

+(PPi) 2 {cQ{m)yu x dsdx = Y^J i . (3.24) 

Jo Jo i=1 

Now we estimate J1-J7 as follows: 

First, by the assumptions (Ai)-(As), or (Ai)'-(yl3), Lemma 3.2, Lemma 3.3, Corollary 
3.4 and Cauchy-Schwarz inequality, we have 

f Ji<C, 

J i^\f o {{cQ(m)) p y dx + C J u 2 dx <C + J jf ((cQ(m))' 9 ) 2 dx, 

'h<C J ((coQimo))^ dx + C J u 2 dx < C, 
Ja < C, 

h < c, 



Jr<c[ [ c {x)c^Q(mf +1 u 2 x dxds + C [ [ {cQ{m)) 2 ^- p - l dxds 
Jo Jo Jo Jo 

Cmaxco(x) f [ c^Q{mf +1 u 2 x dxds + C max^m)) 7 ^ 1 / max(cQ(m)) 7 (is 
io Jo Jo 



< 

< C. 

(3.25) 

Now we estimate Jq, which is divided into two cases: 

Case 1. When the initial masses connect to vacuum continuously (corresponding to the 
boundary condition (2.10)), we have Jq = 0. 

Case 2. When the initial masses connect to vacuum discontinuously (corresponding to 
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the boundary condition (2.8)), we have by Young inequality and Lemma 3.5 



J 6 = -(Ppi) 2 / (cQ(m)y-P(l,s)((cQ(m))P(l,s)u(l,s))ds 



+((3 Pl ) 2 j\cQ(m)y-P(0,s) (cQ(m)f(0,s)u(0,s)) ds 

<C J* {|(cQ(m)) n/3 (l, s)u n (l, s)\ + |(cQ(m)) n/3 (0, s)u n {0, s)\ } ds 

+C J* {(cQ(m)) (7 - /3) ^r(l,s) + (cQ(m)) (7 -^^(0,s)} ds 

<C + C I (cQ(m)) nl3 (-,s)u n (-,s) ds. 
Jo i°°([o,i]) 

Substituting (3.25) and (3.26) into (3.24), we have for Case 1 and Case 2 



(3.26) 



\ I! (ww)0> + wrw I! I! ( (cw) ^)l dsdx 



<C + J 8 , 



(3.27) 



where 



Js = C 



L°°([0,1]) 



(cQ(m))^(;s)u n (;S) 

By the embedding theorem W 1 ' 1 ^, 1]) <^ L°°([0, 1]), we have 

Js<C f [\cQ(m)) nfS u n (x,s)dxds + C f C ((cQ(m)) nfS u n (x, s)) dxds 
Jo Jo Jo Jo v ' x 

{cQ(m)y U 2n (x,s)dxds + C [ [ (cQ(m)) 2n ^dxds 

Jo Jo 

[ [ nP(cQ(m)) n > 5 - 1 (cQ(m)) x u n (x,s)dxds+ [ [ nicQim^u^u^xds 
Jo Jo Jo Jo 

<C J max(cQ(m)) 7 u 2n dx^jds + C [ max(cQ(m)) 7 f (cQ(m)) 2n ^dxds 



< C 



+ 



o [o,i] 

+C f f {cQ{m)) 2n ^- p u 2n dxds + %P * f [ {{cQ{m))^) 2 dsdx 
Jo Jo (P + 7) Jo Jo v Jx 



'o Jo 

rt /•! 



+C [ [ (cQ(m) f +1 u 2n - 2 u 2 x dxds + C [ [ {cQim)) 2 ^^ 1 dxds 
Jo Jo Jo Jo 

<C + max(cQ(m)) 2n/3 - 2 T- /3 / max(cQ(m)) 7 { / u 2n dx ] ds 

Jo I - 1 ] \Jo ) 



1iH 2 pi 



((cQ(m))^) 2 dsdx + Cmax(c (x)) J J J ] Q{mf +1 u 2n - 2 u 2 x dxds 



(/3 + 7) 2 

+Cmax(cQ(m)) 2n/3 " /3 ^ 1 - 7 J max(cQ(m)y ds 



l /•< 



<C + 



(P + l) 2 Jo Jo 



(cQ(m))V) 



dsdx. 



Here we have used Lemma 3.3, Corollary 3.4, Lemma 3.7, and n > 2 ^^ . 
Substituting (3.28) into (3.27), we get (3.21). This proves Lemma 3.8. 



(3.28) 
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3.2 Asymptotic behavior of cQ(m) 

To apply the uniform estimates obtained above to study the asymptotic behavior of the 
mass functions m(x,t), n(x,t) with respect to the time t, we introduce the following lemma. 
The proof is quite simple and the detail is omitted. 



Lemma 3.9. Suppose that g(t) > for t > 0, g(t) G L^O, oo) and g'{t) G L^O, oo). Then 
lim = 0. 

t— >oo 

Now we prove Theorem 2.2. Let 

S(t)= / (cQ(m)) 7 (x,i)(ix. (3.29) 
Jo 



10 

Integrating (3.9) with respect to x over [0, 1], we have 

rt pi 

'o Jo 



[ [ (cQ(m)y(x,s)dxds <C, (3.30) 
Jo Jo 



which implies g(t) G ^ 1 (0, oo). 

Now we prove g'(t) G L 1 (0, oo). By the second equation of (2.7) and using Cauchy-Schwarz 
inequality, we obtain 



POO PC 

/ \g'(t)\dt = 1 
Jo Jo 



1 

(cQ(m)y- l c(x)Q(m) t dy 

o 



dt 



l 

7P/c° ' Q(my +1 u x dx 

o 



dt 



f'OO r\ /*00 rl 

C / / c^Q(m) 1+ ^u 2 x dxdt + C / / c^-^Q(m) w ~^dxdt. (3.31) 
Vo jo jo Jo 



< 



By (3.8), (3.9), and the assumptions (A\) or {A\)' , we can estimate the last term on the 
right-hand side in (3.31) as follows: 

f f 1 c 2 ^Q(m) w -^dxds < f max(cQ(m)y C c 7 "^ 'Q(m) 7+1 ^ ''dxds < C. 
Jo Jo Jo I - 1 ! Jo 

Substituting the above inequality into (3.31) and using Lemma 3.2, we deduce g'{t) G 
oo). 
Consequently, 

lim g(t) = 0. (3.32) 

t^oo 

By (3.32), Lemma 3.3 or Holder inequality, we have 

lim f (cQ(m)) x (x,t)dx = 0, (3.33) 
t^ooJ Q 



for any < A < oo. 

Now we prove Theorem 2.2, which is divided into two cases: 
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Case 1. When the initial masses connect to vacuum continuously (corresponding to the 
boundary condition (2.10)), choosing k > (3 > and applying (3.33), Lemma 3.8 and Holder 
inequality, we have 



X 



< (cQ{m)) k = I {{cQ{m)f ) y dy 

fc(cQ(m)) fc - /3 (cQ(m)) /3 - 1 (cQ(m)) y dy 



o 



o 



x 



| / (rQ(,»)/- ; ((rQ(/»))''M,/ 



o 



1 . \3 



<C( {cQ{m)f k -^dx / ({cQ{m)Yf x dx 



/ vo 

"l \ a 

< C ( / (cQ(m)) 2fc " 2/3 cfe) -»■ 0, as t ->• oo. 



Case 2. When the initial masses connect to vacuum discontinuously (corresponding to the 
boundary condition (2.8)), choosing k > j3 > and applying (3.33), Corollary 3.6, Lemma 
3.8 and Holder inequality, we have 

< (cQ(m)) k = (cQ(m)) k (0,t) + f \{cQ(m)) k ) y dy 



o 



<C(l + tpF*+( jf (cQ(m)) 2k ^d x y i^J {{cQ(m)f) 2 x dx 



i 



<C{l + t)~T=P + Clj {cQ{m)) 2k - 2l3 dxj -+0, ast^oo. 

Combining the above two cases, we have 

(cQ(m))(x, t) — >• 0, as t — > oo, 



which implies 



Thus 



and 



n m ,. n 
hm — • = lim = 0. 

t^oo m pi — m t^oo pi — m 



Tl 

lim n(x, i) = lim • (p/ — m) = 0, 

t->oo i— >oo pi — m 



lim m(x, i) = 0, 



t— >oo 



for any x G [0, 1]. 

This completes the proof of Theorem 2.2. 

4 Decay rates of the mass functions 

Now we are in the position to estimate the stabilization rates of the mass functions 
m(x,t), n(x,t) as t — > oo. 

To do this, introduce a new function w(x,t) defined as follows (cf. |13j): 

w ( x ,t) = PlU (x,t)- T ±- t l Q^dy + j^l I Q^dyd*. (4.1) 
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By direct calculation, we have 

w x = piu x - | (4.2) 
(1 + t)Q(m) 

and 

w 

w t + = Pi u t- ( 4 -3) 

Here we have used the fact that 

/ u(x,t)dx = / uo(x)dx, 
Jo Jo 

(see (3.1)). Assume now uo(y)dy = for the simplicity of presentation. 
Thus the auxiliary functions w and Q(m) satisfy the following 

' d t c = 0, 

d t Q(m) + Q{mfd x w + = 0, (4 . 4) 

Then we have 

Lemma 4.1. Let (c(x), u(x, t), Q(m)(x, t)) be a global weak solution to the fixed boundary 
value problem (2.7), (2.8), (2.9), or (2.7), (2.10), (2.11). Then for any (3 > 0, j > 1 + 0, 
the following estimates hold for any t > 0, 

Case I: < /3 < 1. 

\l + t) e I' w 2 dx+ il + t)9 ' 1 C c^Q(mf- l dx+ Pl{l + t)e C c^Q{my~ l dx 
2 Jo 1 - P Jo 7 - 1 Jo 

+ - ^ + ^ ™ 2 <^s + jf (1 + s) e c p Q{m) l+ ^w 2 x dxds 

+ ^4 !\^ + s) 6 ~ 2 C cPQimf^dxds + p C' 1 ' 6 !\l + s) e - 1 C c<Q{my- 1 dxds 
1 - P Jo Jo 7 - 1 Jo Jo 



<C, 
where 



(4.5) 

9 = min{ 7 - 1, P} = p. (4.6) 



Case II: = 1. 



hl + t) d C w 2 dx+ Pl{l + t)d I' 1 C iQ{my- l dx 
2 Jo 7 - 1 Jo 

+ i 1 ~ ^) Jo ^ + S ^ _1 / u,2(i:rds + J Q + s )° J Q cQ(m) 2 w 2 x dxds 

+Pi 1 ~ 1 ~ 9 f(l + sf- 1 I' JQimy-'dxds 
I- 1 Jo Jo 

<C + Cln(l + t), (4.7) 
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where 6 = 1. 

Case III: > 1. 



2 Jo 2 U ~ 1 ) Jo 

+ i 1 ~ 2~) + 5)91 jf w2(ixds + + s)9 jf ^Q(m) 1+fi wldxd8 



+ Pl{ l~ l ~ 6) fil + s)'- 1 C^Q{m)^dxds 
< C + C{\n{l + t))\ 



where 



2, 

k 7-/3 : 



for ^>2, 
for ^<2, 



(4.8) 



(4.9) 



and Z = 0, when > 2, whereas / = 1, when < 2. 

Proof. Multiplying (4.4)3 by w, integrating the resulting equation with respect to x over 
[0, 1], using integration by parts, we obtain by the boundary conditions (2.8) or (2.10) 



1 d f 1 2j If 1 , 
— r / w 2 dx-\ / "- 

2 dt J 1 + tJo 



w dx 



f (c^Q(mf +l w x ) x wdx + -j- !\^Q(mf) x wdx- Pl [\c^Q( 
Jo t + t Jo Jo 



m)' y ) x wdx 



= c^Q{mf +l w x w 



+ - — - c^Q(m)^w — pic 1 Q{m,yw 



~ f cPQ(m)P +1 wldx - — *— / cPQ{mfw x dx + p ; / c 1 Q(m)~< w x dx 
Jo l + tJo Jo 



= - / c^m^+^dx 
Jo 



— — - / c?Q{mfw x dx + pi [ c y Q(m)' y w x dx, 
1 + 4 Jo Jo 



(4.10) 



i.e., 



1 d /- 1 



-— / u; 

2dt 



dx H 3— / w 2 dx + / cPQ{m)^ +l w 2 x dx 

1 + * Jo Jo 

= — / c l3 Q(m) l3 w x dx + pi / c / Q(m) 1 w x dx. 

1 + t Jo " Jo 



(4.11) 



Now we will prove (4.5), (4.7) and (4.8). 

Case I: < (3 < 1 (The proof of (4.5)). 

Notes that 



1 



w x = piu x - 



(l + t)Q(m) \Q(m)J t (l + t)Q(m)' 
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Thus we can estimate the first and second terms on the right-hand side in (4.11) as following: 



1 f 1 

/ c l3 Q(m) l3 w x dx 

= ~TT~tl c ^ H 1(w) ( "(HljQH 



(i -£)(! + *) 

and 



f c^{Q{mf- l ) t dx + f 1 cPQ{mf- l dx, (4.12) 

Jo \ L + t ) Jo 



Pl [ cTQimr^dx = Pl £ cW { (-Q^) t ~ {1 + t)Q(m)} dx 

pi pi 
= / c 7 (Q(m) 7_1 ) t cte - / c<Q{my- l dx. (4.13) 

1 - 7 Jo 1+tJo 

Substituting (4.12) and (4.13) into (4.11), we get 

d [Wv? _a 7g(m)7 -A dx + J_ f 1 w 2 dx 
dtj \2 7-1 ; i + ty 

Jo 1 + Wo 

= { s-r ){i + t) f + (TTJ? f ^("■^'*- <« 4 > 

Multiplying (4.14) by (1 + t) e for some 9 to be determined later, we deduce for any < (3 < 1 

+ f 1 w 2 dx+ {1 ^ t)e ~ 1 C c^Q{mf~ 1 dx+ Pl{l+ f C^Qimy-'dx 
dt yi J 1 - p Jo 7 - 1 Jo 

+ (l - (1 + i)* -1 jf * w 2 dx + (1 + t) 9 J* cPQ{m) l+ Pw 2 x dx 



t)^ 2 C c l3 Q(mf- 1 dx + pi- — - — -{l + tf- 1 C c<Q{my- l dx 
Jo 7-1 Jo 



= 0. 

(4.15) 

Taking 9 = min {/3,7 — 1} = (3 in (4.15) and integrating (4.15) with respect to t over [0, t], 
we deduce (4.5). 
Consequently, 

f c 1 'Q(my- 1 dx < C(l + t)~ e . (4.16) 
Jo 

Case II: /3 = 1 (The proof of (4.7)). Under this case, the first term on the right-hand 
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side in (4.11) can be rewritten as 

cQ(m)w x dx 



1 n 



l + t 



o 



i r 1 „, . r ( i 



(i + O jo 

Similar to (4.15), we have: 



1+Wo IVQM/t (l + t)QM 

/ c(ln(Q(m)) f dx + 1 / c (a;)dx. (4.17) 
jo U + ^J jo 



+ ( 1 + t ^ e i f w2dx + ( l + t ) e J cQ(m) 2 w 2 x dx 

+pi 1 ~ 1 ~ 6 {i + tf- 1 C c<Q{my- l dx 

7 - 1 Jo 
= j t ^(l+t) e - 1 j\\n(Q(m))dx^+(l + t) e - 2 j\ {x)dx 

+ {l-e)(l + t) 8 - 2 [ c\n{Q(m))dx. (4.18) 

JO 

By using lnx < x — 1 for any x > and Lemma 3.3, we have 

/ clnQ(m)dx < / cQ(m)dx < C 
Jo Jo 

and the assumption (A\) or (A\)' implies that 

/ co(x)da 
Jo 



Jx< C. 



Taking 6 = 1 in (4.18) and integrating (4.18) with respect to t over [0,t], we deduce (4.7). 
Consequently, 

f c r Q(m)"'- 1 dx <C(l + t)- e ln(l + t). (4.19) 
Jo 

Case III: (3 > 1 (The proof of (4.8)). 

Rewrite (4.15) as 

+ ^1 _ ^ (1 + t) 9 - 1 J* w 2 dx + (1 + t) e C PQ(m) l+ Pwldx 

+p l 1 ~ l ~° (l + t) - 1 f 1 cVQimy-'dx 
7 - 1 Jo 

= j t [ {l+ p t } (> l 1 f <^Q(m)^ -1 dx| + + t) e - 2 jT 1 JC&mf^dx. (4.20) 
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Integrating (4.20) with respect to t over [0, t], we have 

-(l + tf [ 1 w 2 dx+ Pl{1 + t)9 C ciQimy^dx 
2 Jo 7 - 1 Jo 

+Pi 1 ~ l ~° (\l + s) e - 1 f 1 c<Q{my- 1 dxds 
7 - 1 Jo Jo 

= 1 f 1 w 2 dx + ^l f 1 c i Q{mQ y-i dx _J_ f 1 c P Q{mo f-i dx 
1 Jo 7 ~~ 1 Jo P ~ 1 Jo 

(l + ^-i f l c p Q{m) ^ ldx+ P z f\ l + s) e-2 f 1 c P Q{m f-i dxds 
l Jo P ~ 1 Jo Jo 



(3-1 

J 



^— C c^Qimof^dx + I f 1 w 2 dx + ^(mo) 7 - 1 ^ + / 1 + / 2 . (4.21) 

— 1 Jo ./o 7 — 1 Jo 



By Young inequality, we have 

-1 rl 



h = + " / c^Q{mf- l dx 
P ~ 1 JO 

If 1 o , (/3-l)9 „ -i (,8-1)0 

/ (cQH^^l + 1) — c(l + 1)^ 1 —^ 
Jo 



/3-Uo 



< pl }} + ^ / ,1 c 7 Q(m) 7 - 1 ^ + C(l + t) ( ^ 1 " I ^ ) ^ C co{x)dx 
2(7 - 1) Jo Jo 

< + ^- c^QM 7 - 1 ^ + C(l + t) {e ~ 1 - Li ^ ) ^, (4.22) 
2 (7 - 1) Jo 



and 



I 2 = ^-4 [\l + s) e - 2 j^c^Qimf-'dxds 



P 1 Jo 

^r 1 :^ f\l + s)^ CciQimT^dx + C l\l + sf'^ ii ^^ds. 
Jo Jo Jo 



2(7 - 1) 

(4.23) 

Substituting (4.22) and (4.23) into (4.21), we have 

til + sf- 1 f w 2 dxds+ [\l + sf f c^Q{m) 1+p wldxds 
o Jo Jo Jo 

fil + sf- 1 C ' c*Q(my- 1 dxds 
2(7 - 1) Jo Jo 

< C + C{l + tf~^ +C I {l + sf'^^ds. (4.24) 

Jo 
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H 1 " 



Taking 6 = 2, when ^ > 2 in (4.24), we have 

f c r Q(my- 1 dx < C(l + ty 6 . (4.25) 
Jo 

Taking 9 = 2=jj, when ^ < 2 in (4.24), we have 

/ C TQ(m)T- 1 dx<C(l + t)- 9 ln(l + t). (4.26) 
Jo 

This completes the proof of Lemma 4.1. 

The proof of Theorem 2.3. Under the boundary condition (2.8), for < f3 < 1 or 

P > 1) ^p^j > 2, choosing some constant 2k = 7 — 1 + 2/3 and using the assumption (Ai), 
Corollary 3.6, Lemma 3.8, (4.16) and (4.25), we have 

(cQ(m)) k (x,t) = (cQ(m)) k (0,t) + f {(cQ(m)) k ) y (y,t)dy 

Jo 

< C(l + t)~T* + C ( y J Q {(cQ(m)f) 2 x dx^j 2 (cQM) 2 ^ 2 ^ 2 



<C(l + t)~^P + ^-c^Q{mf- l dxy 
< C(l + t)~^ +C(l + t)~3 

<C(l + t)-i, (4.27) 

which implies 

(cQ(m))(x,t) < C(l + t)~^ = C(l + t)~^W, (4.28) 

i.e., 

n(x '^ - < C(l + t)~^W. (4.29) 
pi-m(x,t) 

Thus 

n(x, i) = • (p, - m(x, t)) < C(l + t)-T&5!5 , (4.30) 
Pi TTlyX,t) 

and 

m(x,i) = n(x,t) • c(x) _1 < C(l + f)"^, (4.31) 

for any 2; £ [0, 1]. 

Similarly, if (3 = 1 or (3 > 1, < 2, we have 

(cQ(m)) k (x,t) < C(l + i)-Vln(l + <), 

which implies 

n(x,i) < C(l + t)~^W(ln(l + t))^W, (4.32) 

and 

m(x,t) < C(l + t)~^W(ln(l + t))^W, (4.33) 
for any x £ [0, 1]. Here we have used (4.19) and (4.26). 
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Under the boundary condition (2.10), for 0</3<lor ( S>l, > 2, choosing some 
constant 2ki = + 2/3 and using the assumptions (Ai)' , Lemma 3.8, (4.16), (4.25) and 
Holder's inequality, we have 

(cQ(m)) k ^x,t) = [ X ((cQ(m))^) y (y,t)dy 
Jo 

< C ^j\(cQ{m)Yf x dx^ 2 (^j\cQ(m)) 2k ^d: 



<(' [ I — - r ■(.){„,) ■' d.r 
/O C2 

< C ( jf c^Q{my- l dx\ * ( [ —^dx^ 4 



which implies 



i.e., 



o c o( x ) 

<C(1 + <)"*, (4.35) 



(cQ(m))(x,t) < C(l + t) **i = C(l + (4.36) 



n(M) <C(l + t)-™?, (4.37) 



pi - m(x,t) 



Thus 



and 



n(x, t) = n(a, f) - • (p, - m(x, *)) < C(l + t)"^M? , (4.38) 
pi-m(x,t) 



m(x,t) < C(l +t) FW, (4.39) 



for any x G [0, 1]. 

Similarly, if /3 = 1 or /3 > 1, < 2, we have 

(cQ(m)) fel (x,i) < C(l + t)-f(ln(l+t))3, 

then 

n(x,t) < C(l + i)~7-i%< (l n (i + £)) 7-1+4/3 ; (4.40) 

and 

m(x,t) < C(l + i)~7-i+4/3 (i n (i + t)) 7-1+4/3 ) (4.41) 

for any x 6 [0, 1]. Here we have used (4.19) and (4.26). 
The proof of Theorem 2.3 is completed. 
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